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We analyze the observations of the X-ray pulsar LMC X-4 performed by the INTEGRAL 
observatory and the All-Sky Monitor (ASM) of the RXTE observatory over a wide energy 
range. The observed hard X-ray flux from the source is shown to change by more than 
a factor of 50 (from ~ 70 mCrab in the high state to ~ 1.3 mCrab in the low state) on 
the time scale of the accretion-disk precession period, whose mean value for 1996-2004 was 
determined with a high accuracy, P wec = 30.275 ± 0.004 days. In the low state, a flare 
about 10 h in duration was detected from the source; the flux from the source increased by 
more than a factor of 4 during this flare. The shape of the pulsar's broadband spectrum is 
essentially invariable with its intensity; no statistically significant features associated with 
the possible resonance cyclotron absorption line were found in the spectrum of the source. 
© 2005 Pleiades Publishing Inc. 
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INTRODUCTION 



The X-ray pulsar LMC X-4 in the Large Magellanic Cloud (LMC) (the distance to the 
object is d = 50 kpc) is a high-mass binary with a pulsation period of ~ 13.5 s (Kelley 
et al. 1983), in which the compact object is eclipsed by its optical companion, an 08 star 
of the 14th magnitude with a mass of 20 M (Chevalier and Ilovaisky 1977), every ~1.408 
days (Li et al. 1978; Lang et al. 1981; Levine et al. 2000). Levine et al. (2000) provided 
the following orbital parameters for the binary: a x sini = 26.333 ± 0.019 light seconds, the 
eccentricity e < 0.003(2<r), and the epoch of the zero orbital phase T =MJD51110. 86571. 

Lang et al. (1981) pointed to the existence of a superorbital period in the binary, ~ 30.5 
days, within which the intensity of the source changes by a factor of about 60. These authors 
also assumed that this effect is produced by the blockage of direct X-rays by a precessing 
accretion disk tilted with respect to the orbital plane, much as is the case in the binary 
Her X-l (Tananbaum et al. 1972). Based on the measurements of almost the same X- 
ray flux observed from the source during eclipses in its low and high states, Woo et al. 
(1995) concluded that the intrinsic luminosity of the pulsar is constant, which also confirms 
the assumption made by Lang et al. (1981). Having analyzed the RXTE and GINGA 
observations, Paul and Kitamoto (2002) estimated the rate of decrease in the precession 
period of the accretion disk to be P ~ —2 x 10~ 5 s s -1 . 

The flaring activity of LMC X-4 has been widely discussed in the literature (see, e.g., 
Epstein et al. 1977; Skinner et al. 1980; Kelley et al. 1983; Levine et al. 1991). A flare with 
a duration from dozens of seconds to dozens of minutes is detected from the source during 
its high state, on average, once a day. However, there is also evidence (Woo et al. 1995) for 
the presence of bursts in the low state. 

In different periods, the pulsar LMC X-4 demonstrates both a spin-down and a spin-up, 
suggesting that its period is close to its equilibrium value. In this case, according to the 
disk-accretion theory (Ghosh and Lamb 1979), a neutron star with an equilibrium period of 
~ 13.5 s and an X-ray luminosity of ~ 4 x 10 38 erg s' 1 would have a very high magnetic 
moment, ~ 10 315 G cm 3 (Woo et al. 1996). 

The spectrum of the source above several keV is described by a power law with a high- 
energy cutoff. The measured absorption column density toward the source (along with 
the best-fit parameters) does not depend on the phase of the precession period; its value, 
N H ~ 5.5 x 10 20 atoms cm~ 2 , is close to the Galactic column density, suggesting the absence 
of strong internal absorption in the binary (Naik and Paul 2003). Different authors (Levine 
et al. 1991; Mihara 1995; Woo et al. 1996; La Barbera et al. 2001) obtained several estimates 
for the presence of a cyclotron absorption line in the source's spectrum in the energy range 
19 to 100 keV, but as yet there are no reliable measurements of its intensity and, accordingly, 
the surface magnetic field of the neutron star. 

In this paper, we present the results of our timing and spectral analyses for the pulsar 
LMC X-4 performed by using INTEGRAL and RXTE observations over a wide (1-100 keV) 
energy range. Preliminary results were obtained by Lutovinov et al. (2004). 
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OBSERVATIONS 

The INTEGRAL International observatory (Winkler et al. 2003), which carries four 
scientific instruments that allow the emission from astrophysical objects to be studied over 
a wide wavelength range (from optical to hard 7 rays), was placed in orbit by a Russian 
Proton launcher on October 17, 2002 (Eismont et al. 2003). 

In this paper, we analyze the data obtained with the ISGRI detector of the IBIS gamma- 
ray telescope (Ubertini et al. 2003) and the JEM-X X-ray monitor of the INTEGRAL 
observatory (Lund et al. 2003). The ISGRI detector has an effective area of ~ 960 cm 2 
at 50 keV, is effectively sensitive to photons in the energy range 20 to 200 keV (the energy 
resolution is ~7% at 100 keV), and allows the image of the sky within a 29°x29° field of 
view (the full coding zone is 9°x9°) to be reconstructed with a nominal spatial resolution 
of ~ 12 arcmin (the angular size of the mask element). See Lebrun et al. (2003) for a more 
detailed description of the detector. The main elements of the two identical JEM-X modules 
(JEM-X1 and JEM-X2) are coded-mask xenonmethane gas chambers placed at about 3.2 m 
from the plane of the detectors. Each module of the telescope has the following technical 
characteristics: the energy range 3-35 keV, the field of view (the full coding zone) 13.2° (4.8°) 
in diameter, and the detector area 500 cm 2 (Lund et al. 2003). 

About 500 individual pointings were made as part of the studies of the LMC region 
performed by the INTEGRAL observatory from January 2 through January 28, 2003, during 
which the pulsar LMC X-4 was within the field of view of the INTEGRAL instruments. 
The total observing time for the source was ~ 10 6 s; at all pointings, it was in the full 
coding zone. We also used data from the All-Sky Monitor (ASM) of the RXTE observatory 
(http://xte.mit.edu/ASM_lc.html) to analyze the long-period characteristics of the source. 
These data are daily averaged fluxes from the source in the energy range 1.3-12.2 keV. 

We constructed an X-ray image and performed a spectral analysis of the IBIS data using 
the methods described by Revnivtsev et al. (2004) and Lutovinov et al. (2003). An analysis 
of a large set of calibration observations for the Crab Nebula revealed a systematic error of 
~10% in the estimated absolute flux from the source over a wide energy range; the spectral 
shape was reconstructed with an accuracy as high as 3-5%. We took this into account in our 
spectral analysis by adding a systematic error of 5%. All of the errors given in this paper 
are purely statistical. The standard OSA software of version 4.1, which was provided by the 
INTEGRAL Science Data Center (http:/ /isdc.unige.ch), was used for our timing analysis of 
the source on pulsation period time scales and for our analysis of the JEM-X data. 

Figure 1 shows the IBIS 20-60 keV map of the LMC region. This map was constructed 
over the entire period of the source's observations in January 2003. Apart from the pulsar 
LMC X-4 itself, the black hole LMC X-l and the single pulsar PSR 054069 are detected 
in the image at a statistically significant level. No emission is detected from the SN 1987A 
remnant (Shtykovskiy et al. 2005). 
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TIMING ANALYSIS 

The Precession Period 

As was said in the Introduction, there is a period of about 30 days in the binary LMC X-4 
that is most likely related to the eclipse of the neutron star's emitting regions by a precessing 
tilted accretion disk. To determine the exact precession period, we used the ASM/RXTE 
data obtained from January 1996 through July 2004. The best precession period of the 
accretion disk was found by the epoch-folding technique (Leahy et al. 1983) to be 30.275 
days. We used the following method to estimate the error of this measurement: for each 
point on the light curve of the source within a la error in the flux, we randomly chose a point 
that was used as a new flux from the source. In this way, we modeled another light curve for 
which, in turn, the best period was determined. We obtained a total of about 100 periods 
whose spread was considered to be a la confidence interval. Such an analysis yielded the 
following mean precession period of the accretion disk for 1996-2004: P prec = 30.275 ±0.004 
days. 

The hard X-ray (20-50 keV) light curve of the source shown in Fig. 2a has a characteristic 
shape. About 2 days after the beginning of its observations, the pulsar switched to the on 
state; in another 5 days, its flux increased by a factor of about 6-8 and reached a peak (~ 70 
mCrab). The arrows in the figure indicate the switch-on and switch-off times estimated more 
than 20 years ago by Lang et al. (1981); the switch-on time was calculated by using the 
ephemerides derived by these authors; the switch-off time was determined by assuming that 
the high state of the source lasted ~ 60% of the entire precession period (Lang et al. 1981). 
We see that the presumed zero phase (switch-on time) of the superorbital cycle is greatly 
shifted relative to the current INTEGRAL measurements. Note that the same relationship 
between the durations of the high and low states is also obtained from the INTEGRAL 
observations of the source (Fig. 2a). 

Unfortunately, the switch-on time of the object under study was not observed by the 
INTEGRAL observatory directly, and we had to use interpolation to determine the zero 
phase of the precession period. The switch-on time of the pulsar LMC X-4 estimated in this 
way is MJD 52644.5 ± 1.0. Owing to good statistics and a long series of observations, we 
were able to determine the mean precession period with a high accuracy (see above), which 
ensures the determination of the presumed zero phase 20 years ago with an accuracy of ~ 1 
day by assuming the period to be constant. Using our value of the precession period, we 
can calculate the presumed zero phase closest to that measured by Lang et al. (1981). Our 
result is either ~ 18.5 days before or ~ 12 days after this time. Thus, the precession period 
of the accretion disk, on average, must either increase with P ~ 1.3 x 10~ 3 s s^ 1 or decrease 
with P ~ —2 x 10~ 3 s s -1 . Here, we give lower limits on the rate of change in the period, 
since we disregard the change that is a multiple of an integer number of periods. 

On the other hand, the mean rate of change in the precession period can be estimated 
by using our values and those obtained by Lang et al. (1981). This estimation yields a 
value of P = (—2.5 ± 0.8) x 10~ 5 s s -1 , which is much lower than our previous estimate 
obtained by extrapolating the switch-on times and agrees, within the error limits, with the 
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value determined by Paul and Kitamoto (2002). 

The observed discrepancy in the estimated rates of change in the precession period may 
stem from the fact that the precession period of the accretion disk for the pulsar LMC X-4 
has evolved not uniformly, but with spinups/spindowns much higher than the average level. 

Using long-term ASM/RXTE observations (January 1996-August 2002), Klarkson et al. 
(2003) constructed a dynamic power spectrum for the source under study. Throughout the 
period of observations, the precession period was constant near 30.28 days. However, due to 
the 3cr error of ±0.46 days, no firm conclusion about the precession-period variations can be 
reached. 

Figure 2b shows the averaged light curve of the pulsar LMC X-4 constructed from 
ASM/RXTE data over ~ 8.5 years of its observations and folded with the best period 
of 30.275 days. Interestingly, there is a statistically significant rise in the source's intensity 
with a duration of about three days near phase ~ 0.8 of this cycle. A similar feature in the 
light curve of the pulsar can be seen in Figs. 1 and 3 from Lang et al. (1981). A similar rise 
in the source's intensity during the off state is also observed in another X-ray pulsar, Her 
X-l (Jones and Forman 1976). 

The Orbital Period 

Figure 3a shows the 20-50 keV light curve for the pulsar LMC X-4 when the source was 
near the peak of its observed X-ray flux. The times of the eclipses of the X-ray source by its 
optical companion are clearly seen from this figure. The ingress and egress were accompanied 
by smooth variations in the pulsar's observed intensity. However, a significant rise in the 
flux from the source with a duration of about three hours was observed a short time (of 
the order of an hour) after the egress. The ingress and egress duration can be estimated 
from the source's light curve, ~ 3 ks. At the neutron star's orbital velocity of about 500 km 
s _1 (Kelley et al. 1983), this corresponds to the transition region's size of ~ 1.5 x 10 6 km, 
which is comparable to the sizes of the atmospheres of late-0 main-sequence stars (Woo et 
al. 1996). 

The time dependence of the hardness, which is defined as the ratio of the fluxes from 
the source in the 40-60 and 20-40 keV energy bands, is shown in Fig. 3b. We see a 
tendency for the radiation from the object under study to soften as it approaches an eclipse, 
but this dependence breaks down immediately after the eclipse. The hardness during the 
eclipses themselves is not shown in the figure, because no statistically significant flux can be 
recorded from the source. 

The dashed lines indicate the eclipse ephemerides taken from the papers by Li et al. 
(1978), Lang et al. (1981), and Levine et al. (2000). Note that the orbital period in the 
binary is highly stable and has been almost constant over more than twenty years of the 
pulsar's studies. The INTEGRAL observations of the source in 2003 confirm this conclusion 
(see Fig. 3). Several authors (Woo et al. 1996; Levine et al. 2000) have attempted to measure 
the rate of change in the orbital period, but the significance of the values obtained is low 
{Porb/Porb = (-5.3 ± 2.7) x 10~ 7 yr" 1 and P orb /P orb = (1.1 ± 0.8) x 10~ 6 yr" 1 , respectively) 
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and is consistent with the observational data. 

Note that the low-state flux from the source does not drop to zero and is recorded by 
the INTEGRAL instruments, although it decreases. Thus, for example, according to the 
observational data spanning the time interval MJD 52662.5-52665, the 20-50 keV flux from 
the source was 2.6 ±0.4 mCrab, while several days later, MJD 52667-52668, it decreased by 
almost a factor of 2, 1.3 ± 0.5 mCrab. 

The intensity of the pulsar LMC X-4 decreases sharply during its X-ray eclipses. We 
analyzed the flux during eclipses for two time intervals, when the source was in its high and 
low states (MJD 52648.8 and MJD 52665, respectively). In both cases, no emission was 
detected from the object under study at a statistically significant level, and we obtained 
only upper limits on its 20-50 keV flux, which proved to be identical, ~ 1.2 mCrab (lcr). 
Note that during the X-ray eclipse on January 20, 2003, (MJD 52659.94-52660.17), when 
the source's intensity decreased, its observed 20-50 keV flux was significantly higher than 
the above upper limits, 3.7 ±1.0 mCrab, but the significance of this measurement is low. 
A similar analysis for the soft (0.1-2.4 keV) X-ray emission was performed by Woo et al. 
(1995). The results obtained by these authors suggest that the source's intensity is constant 
during its eclipses, which they interpreted as evidence that the pulsar's intrinsic luminosity 
is constant. This conclusion was based on the assumption that the X rays observed during 
an eclipse are scattered by the widely distributed matter around the source, possibly the 
corona above the outer edges of the accretion disk, mush as is the case in the pulsar Her X-l 
(Lutovinov et al. 2000). 

Flares 

Many authors have described the flaring activity of the pulsar LMC X-4. Short flares (up 
to 1000 s in duration), during which the pulsar's intensity changes by a factor of 2 to 5, are 
most characteristic of the source under study (Levine et al. 2000). 

An rise in the observed flux from the source was recorded during its low state (MJD 52666) 
on a time scale atypical of it. The flare duration was about 10 h (orbital phase 0.3-0.6), and 
the peak intensity (~ 10 mCrab) exceeded the average low-state level by a factor of about 4. 
Figure 4 shows the flare profile constructed from IBIS data in the 20-50 keV energy band. 

SPECTRAL ANALYSIS 

In the previous section, we showed that a significant number of features could be distin- 
guished in the light curve of the pulsar LMC X-4. Therefore, studying the spectral evolution 
of the source with time is of considerable interest. Another goal of the spectral analysis is the 
search for the possible cyclotron resonance absorption line whose existence was mentioned 
by several authors (see the Introduction). 

A power law with an exponential high-energy cut-off (White et al. 1983), which is typical 
of this class of objects, was chosen as a basis for the spectral analysis of the source over a 
wide (4-100 keV) energy range: 
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where E is the photon energy in keV, A is the normalization of the power-law component, 
a is the photon spectral index, E c is the cutoff energy, and Ef is the e-folding energy in the 
source's spectrum. We identified a total of four states in which an independent analysis of 
the pulsar's spectrum was performed: the high state (without including the X-ray eclipse 
times), the bursts immediately after the source's egress, the low state (without including the 
X-ray eclipse times), and the flare occurred on January 27, 2003 (MJD 52666). 

The best-fit parameters for the source's spectra based on model (1) are given in the table. 
We were able to reconstruct the spectra over a wide energy range (using JEM-X data) only 
in the high state. Since the JEM-X monitor is not sensitive enough, the source was detected 
at a statistically significant level only by the IBIS telescope. Therefore, when fitting the 
spectra in this state, we fixed the parameters in the standard X-ray energy range (< 20 keV) 
at the values obtained for the high state. We see from the table that the shape of the source's 
spectrum and its parameters (slope and cutoff parameters) remain almost constant with the 
pulsar's intensity variations. Figure 5 shows the energy spectrum of LMC X-4 reconstructed 
from INTEGRAL data in a wide energy range for the high state. 

To test the hypothesis about the presence of a cyclotron feature in the source's spectrum 
within the energy range 4-100 keV, we modified the fitting model by adding the correspond- 
ing component: 



exp 



A cyc W* yc {E/E cyc y ~ 
\E-E cyc f + Wl 



eye 



(2) 



where E is the photon energy in keV, E cyc is the cyclotron energy, W cyc is the width of 
the cyclotron line, and A cyc is its depth. 

Using the modified model, we fitted the pulsar's spectrum during the high state, except 
the X-ray eclipse times. We used the same procedure that was employed to study the X- 
ray pulsar KS 1947+300 (Tsygankov and Lutovinov 2005): the energy of the center of the 
presumed cyclotron line E cyc was varied over the range 5-100 keV at 5-keV step, while its 
width was fixed at 5 keV; for each trial energy of the line center, we tested the significance 
of an improvement in statistics using the A% 2 test. As a result, we found no such energy of 
the cyclotron line E cyc in the energy range under study whose inclusion in the model would 
lead to an improvement in the quality of the spectral fit by more than 2a. 

The result obtained may suggest that either the energy of the cyclotron line lies outside 
our energy range (4-100 keV) or the INTEGRAL instruments are not sensitive enough for 
the cyclotron line to be detected in the spectrum of the source LMC X-4. The relationship 
between the spin period of the neutron star and its luminosity, assuming that the period is 
close to its equilibrium value, argues for the hypothesis of a strong magnetic field (> 10 13 
G) in the binary (Woo et al. 1996). 



-7- 



CONCLUSIONS 



We presented the results of our long-term timing and spectral analyses for the X-ray 
pulsar LMC X-4 performed by using the INTEGRAL observations in January 2003 and 
the long-term (1996-2004) ASM/RXTE observations. We determined the precession period 
(30.275 ± 0.004 days) averaged over the last ~ 8.5 years with a high accuracy. Since the 
source's switch-on times predicted by different authors do not coincide, we showed that this 
parameter is not constant and most likely varies nonuniformly. The IBIS/INTEGRAL data 
revealed small bursts that emerged after the egress of the X-ray source. We detected a 
flare about 10 h in duration during the low state. On the time scales of the orbital period, 
its value determined more than 20 years ago is stable and satisfies well our observational 
data. We performed a spectral analysis for different states distinguishable by the object's 
intensity. The spectrum of the source is described by a power law with a high-energy cutoff, 
which is characteristic of this class of objects. To test the hypothesis about the presence of 
a cyclotron feature in the source's spectrum within the energy range 4-100 keV, we properly 
modified the fitting model. This analysis showed that there is no such feature in the source's 
spectrum at a confidence level exceeding 2a. A comparison of the derived constraints with 
the results of other authors (Woo et al. 1996) is more likely indicative of a strong magnetic 
field (> 10 13 G) on the neutron-star surface than a weak field (< 5 x 10 11 G). 
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Best-fit parameters for the spectrum of LMC X-4 derived from INTEGRAL data (JEM-X + IBIS) 



State 


Luminosity, 


a 


E c , 


Ef, 


X 2 n(N) c 




y 1 n 37 pra h 

aiu ei s 




Ke v 


Kt; V 




High (MJD 52646-52650) 


38.7 


0.20 ±0.15 


9.1 ±0.8 


11.0 ±0.6 


0.88(128) 


Flares (in high state) 


47.9 


0.44 ±0.20 


8.9 ±1.0 


11.9 ±0.8 


0.82(126) 


Low (MJD 52662.5-52668) 


1.4 


0.2 d 


9.1 d 


14.8 ±3.2 


1.19(4) 


Flare (MJD 52666) 


3.1 


0.2 d 


9.1 d 


17.1 ±5.3 


0.41(6) 



All errors are given at a la level. 

The 4-100 keV luminosity at an assumed distance to the source of 50 kpc. 
The %2 value normalized to the number of degrees of freedom N. 
The parameters were fixed at the values obtained for the high state. 
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Fig. 1: IBIS/INTEGRAL 20 60 keV image of the LMC region containing the X-ray 
pulsar LMC X-4. The sources detected at a statistically significant level are labeled; 
the position of the SN 1987A remnant is indicated by the cross. 
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Fig. 2: (a) Hard X-ray (20 50 keV) light curve for the pulsar LMC X-4 constructed 
from IBIS/INTEGRAL data. The arrows indicate the switch-on and switch-off times 
for the source, as estimated by other authors, (b) The light curve for the source 
constructed from ASM/RXTE data over eight years of observations and folded with 
the best period derived from these data. The errors correspond to one standard 
deviation. 
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Fig. 3: (a) Orbital light curve for the pulsar LMC X-4 constructed from IBIS data 
during its high state. The dashed lines indicate the ingress and egress times estimated 
by Lang et al. (1981). (b) Variations of the source's hardness on the time scale of the 
orbital period. The errors correspond to one standard deviation. 
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Fig. 4: Flare detected by the IBIS telescope of the INTEGRAL observatory from the 
pulsar LMC X-4 on January 27, 2003, during the low state. The solid line indicates a 
Gaussian fit to the flare profile. 
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Fig. 5: Energy spectrum of the source LMC X-4 during its high state constructed 
from JEM-X and IBIS data. The dots indicate the experimental spectrum; the solid 
line represents its power-law fit with a high-energy cutoff. 



